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Abstract 

We present the heavy flavour mass corrections to the order a s corrected longitudinal {(t£) 
and transverse (<tt) cross sections in e + e~-collisions. Its effect on the value of the running 
coupling constant extracted from the longitudinal cross section is investigated. Furthermore 
we make a comparison between the size of these mass corrections and the magnitude of the 
order a 2 s contribution to <jt and ul which has been recently calculated for massless quarks. 
Also studied will be the changes in the above quantities when the fixed pole mass scheme 
is replaced by the running mass approach. 
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Experiments carried out at electron positron colliders like LEP and LSD have provided us 
with a wealth of information about the constants Jl| appearing in the standard model of the 
electroweak and strong interactions. One of them is the strong coupling constant a s which can 
be extracted from various quantities. Examples are the hadronic width of the Z-boson, the total 
hadronic cross section, event shapes of jet distributions and jet rates (for a review on this subject 
see 0). Another quantity from which a s can be extracted is the longitudinal cross section 
which is measured in the semi inclusive reaction 

e~ + e + -> V -> P + "X" . (1) 

Here V represents the intermediate vector bosons Z and 7 and "X" denotes any inclusive final 
state. Furthermore P stands for a hadron which is produced via fragmentation by a quark or a 
gluon. As we will see below the perturbation series in QCD for <jl starts in order a s provided 
the quark which couples to the vector boson is massless. Therefore this quantity will become 
very sensitive to the value of the strong coupling constant. Some time ago the perturbation 
series was only known up to order a s (see 0, However it turned out that the result for 

the longitudinal cross section was much smaller than its experimental value measured at LEP 
H [I| which indicated that the higher order corrections may become rather large. The latter was 
confirmed by the second order corrections, computed very recently in [^], which amount to about 
35 % with respect to the lowest order contribution to 0£. If the second order corrections are 
included the longitudinal cross section agrees now very well with experiment so that it can be 
used to extract a s (see ||). Until now <jl was only evaluated for massless quarks. This will be a 
correct approximation for the light quarks like u, d, s and probably also for c but it will certainly 
fail for the 6-quark, leave alone for the top quark. Therefore one has to compute the mass 
corrections to this quantity at least up to first order in a s . This contribution will be presented 
in this paper and we discuss its effect on the value of a s when extracted from experiment via a^. 

The cross section corresponding to process (|1|) where the hadron P emerges directly, or 
indirectly via the gluon, from the heavy quark H is given by 



ciV H - p (x,Q 2 ) 3, 2/A ^ p (z,Q 2 ) , 3 , 2n da*' p (x,Q^ 
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dx d cos u 8 v J dx 4 dx 
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From the equation above one can extract the transverse (<tt), the longitudinal (<jl) and the 
asymmetric (<ta) cross sections. Further 9 denotes the angle between the outgoing hadron P and 
the incoming electron. The Bj0rken scaling variable is defined by 



^ = 1^ g 2 = Q 2 >0, 0<a;<l, (3) 



where the momenta p and q correspond with the outgoing hadron P and the virtual vector 
boson (Z, 7) respectively. To obtain the total transverse and longitudinal cross sections one has 
to multiply Eq. (0) with x. After integration over x and summation over all hadron species P 
emerging from the heavy flavour H one obtains 
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The result above can be written as follows 



vf(Q 2 )=Wv(Q 2 )K(p) + a AA (Q 2 )h a k (p) with p 



Am 

'Q 2 



(5) 



where m denotes the heavy quark mass. The above definitions for the total cross sections should 
not be confused with the ones which are obtained by integrating d a^jd x without muliplication 
by x. The heavy flavour contributions to the latter quantities were computed for the first time 
in U (see also ]T^|-[r2]). The difference between these two definitions will be pointed out below. 
Finally the total cross section for heavy flavour production is given by 



In Eq. 
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the pointlike cross sections are defined by 
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where N denotes the number of colours corresponding to the gauge group SU(N) (in the case 
of QCD one has N = 3). Furthermore we adopt for the Z-propagator the energy independent 
width approximation 
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In Eqs. (0), (|J) the charges of the lepton and the quark are denoted by 
and the electroweak constants are given by 
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and e„ respectively 



r 1 

Ca,u = —Caa = —Caa 

8 
3 



-C 



Al 



4 sin 2 6 W ) 



(10) 



c 



V,u 



c 



Al 



sin 2 6 



W ) 



V,d 



-C 



Al 



sin 2 ^ 



The functions h l k (k = T, L, I — v, a) in Eq. @ can be obtained order by order in perturbative 
QCD from the singlet quark and the gluon coefficient functions in the following way 



dx x C l k ,s Jx, p, Q 2 / p 2 ) + - 



1-/3 



dxxC l k s Jx,p,Q 2 /p 2 



Here p stands for the factorization as well as the renormalization scale. The result in Eq. (O) has 



been derived from the differential cross section d c^' P / d x (see e.g. 0). The latter can be written 
as a convolution of the fragmentation densities Df and the coefficient functions C^j (i = q, g) (see 
Eq. (2.4) in [0). Because of the momentum sum rule satisfied by the fragmentation densities 
(see Eq. (2.9) in 0) expression (|Tl"D follows immediately. Notice that the functions h l k differ 
from the functions fj. computed in [JT2J (see also the functions Hi (i = 2,6) in JTOj). The latter 
were derived from the first moment of the non-singlet quark coefficient function and they differ 
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from the former which are given by the second moment of the singlet quark and gluon coefficient 
functions presented in Eq. ([TTJ) . However because of Eq. (^j) the sum of the transverse and 
longitudinal components of both functions have to lead to the same perturbation series of the 
total cross section so that one has the relation h T + h l L = f l T + f l L (I = v,a). The functions h l k 
can be expanded in the strong coupling constant a s as follows 



h l k ( P ) = E 

n=0 
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h l k (n) (p). 



The lowest order contributions corresponding to the Born reaction 

V^H + H, 



(12) 
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with V = 7, Z are given by 
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where p is defined in Eq. (|5|). The next-to- leading order (NLO) contributions originate from the 
one-loop virtual corrections to reaction (|13D and the gluon bremsstrahlungs process 



V -^H + H + # 



(15) 



We have computed the order a s contributions to the quark C£ and gluon C k ' g coefficients func- 
tions for the case m ^ in Eq. (fTID and confirmed the result obtained in appendix A of ||. 
Notice that the quark coefficient functions are derived from the process where the hadron P 
emerges from the quark H. The gluon coefficient function corresponds with the process where 
the hadron is emitted by the gluon. After performing the integral in Eq. (O) we obtain the 
following results for k = T,L 
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with 

t = WW (20) 

and the colour factor Cp is given by Cf = {N 2 — 1)/2N. The functions Fi(t) appearing above 
are defined by 

F 1 (t) = Li 2 (t 3 )+4C(2) + iln 2 (t) + 31n(t)ln(l + t + t 2 ) (21) 
F 2 (t) = Li 2 (-t 3 / 2 ) - Li 2 (t 3 / 2 ) + U 2 {-t 1 ' 2 ) - U^t 1 ' 2 ) + 3C(2) + 2 hx{t) ln(l + y/t) 

-2 ln(t) ln(l - Vi) + - ln(t) ln(l + t - VI) - - ln(t) ln(l + t + V~t) (22) 

F 3 (t) = Li 2 (-t) + ln(t)ln(l-t) (23) 
F 4 (t) = 6m(t) - 81n(l -t) -41n(l + t) , (24) 

where C( n )> which occurs in the formulae of this paper for n = 2,3, represents the Riemann (- 
function and Li 2 (x) denotes the dilogarithm. Using Eqs. flTED - QTP| ) one can check that the order 
a x contribution to the total cross section in (||) is in agreement with the literature [|HJ. The 



next-to-next-to-leading order (NNLO) contributions come from the following processes. First 
one has to calculate the two-loop vertex corrections to the Born process ( |i~3"D and the one-loop 
corrections to ([15]). Second one has to add the radiative corrections due to the following reactions 

V ^E + E + g + g (25) 
V^H + H + H + H (26) 

V ^E + E + q + q (27) 

/ (2) 

where q denotes the light quark. The results for h£ presented below are only computed for 
those contributions containing Feynman graphs where the vector boson V is always coupled to 
the heavy quark H so that the light quarks can be only produced via fermion pair production 
emerging from a gluon splitting. Because we are only interested in the ratios 

R *(Q 2 ) = 4lS-, for k = T > L ( 28 ) 

we do not consider other contributions which drop out in the expression above provided we put 
the mass of H equal to zero in the second order correction. The contributions which can be 
omitted are given by one- and two-loop vertex corrections which contain the triangular quark- 
loop graphs |I5| , fll6"| . They only show up if the quarks are massive and are coupled to the 



Z-boson via the axial-vector vertex. Notice that one has to sum over all members of one quark 
family in order to cancel the anomaly. Further we exclude all contributions from reaction ( [ZTD 
which involve Feynman graphs where the light quarks q are coupled to the vector boson V. 
Notice that interference terms of the latter with diagrams where heavy quark are attached to 
the vector boson vanish if the heavy quark is taken to be massless provided one sums over all 
members in one family. 

The order a 2 contributions to the quark and gluon coefficient functions have been calculated 
in 0, [14]. Because of the complexity of the calculation of these functions the heavy quark mass 
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was taken to be zero. This approximation is good for the charm and bottom quark but not for 
the top quark as we will see below. Substituting these coefficient functions in the integrand of 
Eq. flllD (here p = 0) we obtain 
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Where the colour factors are given by Ca = N and Tf = 1/2 (for Cp see below Eq. (|20"D ). 
Further n/ denotes the number of light flavours which originate from process (27). Finaly \i 
appearing in the strong coupling constant a s and the logarithms in Eq. ([30"D represents the 
renormalization scale. Notice that the coefficient of the logarithm is proportional to the lowest 
order coefficient of the /3-function. The logarithm does not appear in hip ^ because h l T = for 



and I 



Since m = there is no distinction between h 
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and for k 



T, L 



anymore unlike in the case of the first order corrections in Eqs. (|iT)|)-(|T9"|) where the heavy quark 
was taken to be massive. Furthermore one can check that substitution of Eqs. (p9|), (|30|) into Eq. 
(|5|) provides us with the order a 2 contribution to the total cross section which is in agreement 
with the results obtained in [ I7| . When the quark is massless we get from Eqs. (|5]), @ the 
following perturbation series for the ratio in Eq. (p8|) up to order a 2 
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We will now show how the results for the above ratios will be modified when the Born and the 
first order contribution are computed for massive quarks. Further we discuss the consequences of 
our findings for the extraction of a s . In order to make the comparison between the massless and 
massive approach to Eq. (]2lf ) we have chosen the following parameters (see [TE\)- The electroweak 
constants are: M z = 91.187GeV/c 2 , T z = 2.490 GeV/c 2 and sm 2 9 W = 0.23116. For the strong 
parameters we choose : A^ IS = 237 MeV/c (jif = 5) which implies a s (Mz) = 0.119 (two-loop 
corrected running coupling constant). Further we take for the renormalization scale fi = Q unless 
mentioned otherwise. Notice that we study for H = c, b at the CM energy Q = Mz- For the 
heavy flavour masses the following values are adopted : m c = 1.5 GeV/c 2 , m;, = 4.50 GeV/c 2 
and m t = 173.8 GeV/c 2 . The results for the bottom quark can be found in table [1|. A comparison 
between the first and second column reveals that on the Born level the diference between the 
massive and massless approach to R T is very small (about two promille). For R\ it is more 
conspicious but the correction due to mass effects, which equals 0.0014 for R L , is still much 
smaller than the order a 2 contribution which amounts to 0.0121. The corrections due to mass 
effects are also smaller than the changes caused by a different choice of the renormalization scale. 



If we choose /i = Q/2or/i = 2<5 one gets R { 1 ] = 0.0509 and R^' = 0.0461 respectively which 
differ by 0.0024 from the central value R L = 0.0485. We also studied the effect of the running 
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m b = 0.0 GeV/c 2 


m b = 4.50 GeV/c 2 


m b {M z ) = 2.80 GeV/c 2 


d(0) 
JrCrp 


1.0 


0.999 


0.999 


JXrp 


0.962 


0.964 


0.963 


d(2) 

JXrp 


0.950 


0.952 


0.951 


R b L 




0.0 


0.0016 


0.0006 


R? 


0.0378 


0.0364 


0.0374 


R? 


0.0499 


0.0485 


0.0495 



Table 1: The ratio Rk = crk/&tot (k = T, L) for 55-pro duct ion. 



quark mass on the ratio in Eq. (|28|) . For this purpose one has to change the on-mass shell scheme 
used in Eqs. (p~6|) - (P~9D into the MS-scheme. This can be done by substituting in all expressions 
the fixed pole mass m by the running mass rn(fi). Moreover one has to add to the first order 
contributions (|1^)-(|19|) the finite counter term 



A^' (1) = m(fj)C F 



4 - 3 In 



m 2 (/i) 



/< 2 



d m 



(33) 

m=m(/i) 



where /i stands for the mass renormalization scale for which we choose fj = Q. Further we adopt 
the two-loop corrected running mass with the initial condition fn(^o) = fi^. Using the relation 
between the MS'-mass and the fixed pole mass, as is indicated by the first factor on the righthand 
side in Eq. (^), we have taken for bottom production /i = 4.10 GeV/c 2 which corresponds 
with a pole mass m b = 4.50 GeV/c 2 . This choice leads to m b [Mz) — 2.80 GeV/c 2 which is 5 % 
above the experimental value 2.67 GeV/c 2 measured at LEP |T9[. The results are presented in 
the third column of table |1|. Comparing the second with the third column we observe that the 
mass corrections decrease because the running mass is smaller than the fixed pole mass. It is 
now interesting to see how the mass terms contributing to R b L affect the extraction of the running 
coupling constant at /i = Q = Mz- To that order we equate the mass corrected formula for R b L 
with the massless result presented in Eq. (|3lD . The latter yields R L = 0.0499 for a s = 0.119 
(see first column, third row of table p]). Using the same number for the massive expression of R^ 
we obtain a s = 0.122 which amounts to an enhancement of 2.5 % for m b = 4.50 GeV/c 2 . In the 
case of a running mass i.e. m b [Mz) = 2.80 GeV/c 2 we get a s = 0.120 so that the enhancement 
becomes 1.0 %. However as we have said before these mass effects are smaller than those caused 
by a variation in the renormalization scale. Following the same procedure we equate Rl at 
different values of /i. Choosing /i = Q one infers from table | that R^ = 0.0499 for a s = 0.119. 

(2) 

This value of R L can also be obtained at /i = Q/2 but then one has to take a s = 0.127. If we 
choose n — 2 Q the result is a s = 0.112. Therefore the uncertainty is about 6% which is larger 
than the mass correction. The latter becomes even smaller when we also add the part coming 
from the bottom quark to the longitudinal cross section consisting of the contributions of the 
light quarks and the charm quark. Concerning the latter quark we want to remark that its mass 
is so small that the massive approach will become indistinguishable from the massless results. 
Since the mass effects up to the order a s level are rather small even for the bottom quark we 
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Rp 




mt = GeV/c 2 


m t = 173.8 GeV/c 2 


m t (Q) = 153.5 GeV/c 2 


d(0) 
JrCrp 


1.0 


0.826 


0.862 


d(1) 

JXp 


0.970 


0.828 


0.834 


d(2) 

JXrp 


0.962 


0.821 


0.826 


Ri 




0.0 


0.175 


0.138 




0.030 


0.172 


0.167 


4 2) 


0.038 


0.180 


0.174 



Table 2: The ratio Rk = Ok/vtot (k = T, L) for tt-production at Q = 500 GeV/c. 

can assume that the second order corrections, derived for m = 0, also apply for m ^ at least 
as long as Q > m. This will be correct at large collider energies for all quarks except for the 
top quark as we will see below. 

For ti production we will study the mass effects up to order a s at a CM energy Q = 
500 GeV/c. Here we have chosen <iq = 166.1 GeV/c 2 which corresponds with a fixed pole 
mass of m t = 173.8 GeV/c 2 . This choice leads to m t (Q) = 153.5 GeV/c 2 . From tabled we infer 
that for this process the mass corrections are huge and they are much larger than the first and 
second order corrections. Therefore the numbers in the last row, presented for R T and R l L , are 
unreliable because they only hold if the mass corrections in the order a 2 s contributions can be 
neglected as was done in Eqs. fl29"|), fl3"0"|). This implies that for the top quark R^ and have 
to be computed for m ^ which will be an enormous enterprise. We also studied the effect 
of the running mass presented in the third column of table |2|. Here the difference between the 
fixed pole mass and the running mass approach is larger than the one observed for the bottom 
quark in table |I|. Another feature is that the order a s correction increases when the running 
mass scheme is used and its sign is reversed with respect to the correction obtained for the fixed 
pole mass approach. Notice that a study of the change in R l k (k = T, L) under variation of 
the renormalization scale makes no sense because of the missing exact result in second order for 
m/0. 

Summarizing our findings we have computed the mass corrections up order a s for the trans- 
verse and longitudinal cross sections which are due to heavy flavours. In the case of charm there 
is no observable difference between the massless and massive approach. For the bottom quark we 
see a difference which leads to an enhancement of the extracted value of a s which is of the order 
of 2.5 % (fixed pole mass m b = 4.50 GeV/c 2 ). or 1.0 % (running mass Mb(M z ) = 2.80 GeV/c 2 ). 
These numbers become much smaller when the light flavour contributions are added to the cross 
sections. A variation in the renormalization scale introduces larger effects on the value for a s 
which are about 6.0 % (see e.g also 0). In the case of top-quark production the mass terms 
cannot be neglected anymore since they are much larger than the second order corrections com- 
puted for massless quarks. 
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